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Abstract 

This  paper  describes  results  of  wind  tunnel  experiments  and  Navier-Stokes  computations  conducted  to  help  understand 
hieh-lift  aerodynamics  about  thin,  fighter  type,  multi-element  airfoils  at  flight  Reynolds  numbers  and  lift  coefficients  up 
to  c,  The  wind  tunnel  tests  were  conducted  in  the  NASA  langley  Research  Center  (LaRC)  Low  Turbulence  Pressure 
Tunnel  (LTPT)  as  part  of  a  cooperative  effort  between  the  Navy,  Boeing,  and  NASA  LaRC.  Surface  pressures,  forces  and 
moments,  transition  data  using  hot  films  were  measured  on  a  two  dimensional  (2-D)  airfoil  model  of  a  Boeing  advanced 
fighter  wing  section  configured  with  a  deflected  leading  edge  flap,  shroud  and  a  slotted  trailing  edge  flap.  Effects  of  Rey¬ 
nolds  number,  trailing  edge  flap  gap  and  overhang,  Gurney  flap,  vortex  generators  on  clmax  were  investigated.  Navier- 
Stokes  computations  were  performed  using  structured/chimera  grid  for  several  of  these  configurations.  The  computa¬ 
tions  were  done  with  the  Baldwin-Barth,  Spalart-Allmaras  one  equation  turbulence  models  and  the  Shear  stress  trans¬ 
port  two-equation  turbulence  model  to  predict  the  Reynolds  number  effects.  The  lift  coefficient  was  predicted  within  3  k 
of  the  experiment  up  to  and  including  stall  angles-of-attack.  Trends  in  maximum  lift  as  a  function  of  Reynolds  number 
and  Mach  number  are  accurately  predicted  indicating  that  the  CFD  method  could  be  used  to  extrapolate  sub-scale  high- 
lift  system  aerodynamic  performance  to  flight  scale  at  angles  of  attack  up  to  and  including  maximum  lift.  The  trends  of 
the  effects  of  flow  control  mechanisms  such  as  Gurney  flaps  were  correctly  predicted.  Results  from  fully  turbulent 
Navier-Stokes  calculations  are  also  correlated  with  the  boundary  layer  transition  data  obtained  from  hot  film  measure- 
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Nomenclature 
Clean  Airfoil  Chord 
Lift  Coefficient 
Maximum  Lift  Coefficient 
Pressure  Coefficient 
Mach  Number 
Ovehang  (%  Chord) 

Reynolds  Number 
Nondimensional  Chord  (%) 
Nondimensional  Span  (%) 

Angle  of  Attack  (Deg.) 

Trailing  Edge  Flap  Deflection  Angle  (Deg.) 
Leading  Edge  Flap  Deflection  Angle  (Deg.) 
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Pitch 

Cpt2 

Cps2 
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Cd 


Shroud  Deflection  Angle  (Deg.) 

Vertical  position  of  the  probe  in  the  test  section 
(in.) 

Model  wake  flow  angle  (Deg.) 

Probe  total  pressure  coefficient 
Probe  static  pressure  coefficient 
Drag  coefficient  at  each  vertical  location 
(zpos2)  in  the  wake  rake  sweep 
Section  profile  drag  coefficient 


Introduction 

The  Navy  depends  on  low-speed  high-lift  aerodynam¬ 
ics,  to  enable  high  performance  multi-role  strike/fight¬ 
er  aircraft  to  operate  from  a  carrier  deck.  Carrier 
suitability  is  always  the  most  challenging  requirement 
which  drives  aircraft  design  and  distinguishes  the  re¬ 
quirements  of  Naval  aviation.  Over  the  last  several  de¬ 
cades  the  approach  speed  of  Navy  high  performance 
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aircraft  has  increased  because  of  additional  multi-mis¬ 
sion  operational  demands  placed  on  the  aircraft.  High 
approach  speeds  have  been  shown  to  be  directly  corre¬ 
lated  to  reduced  operational  safety  and  increased 
maintenance  and  life  cycle  cost.  So,  it  has  become 
very  important  to  study  and  improve  the  high-lift  sys¬ 
tems  in  order  to  reduce  the  approach  speeds  and  im¬ 
prove  safety  of  the  flight  operations  operating  from  a 
carrier  deck. 


The  Recent  advances  and  implementation  of  advanced 
turbulence  models  have  prompted  validation  studies  of 
Computational  fluid  dynamic  (CFD)  codes  of  multi¬ 
element  high-lift  airfoil  flows.  All  of  this  effort  has, 
however,  been  focused  on  the  validation  of  CFD  codes 
and  performance  of  turbulence  models  on  commercial 
type  configurations.  There  are  two  major  geometric 
differences  that  distinguish  modem  high  performance 
multi-role  strike/fighter  military  airfoils  from  com¬ 
mercial  configurations:  1)  leading  edge  shape,  and  2) 
airfoil  thickness.  Integration  of  stealth  requirements 
typically  dictates  sharp  leading  edges  and  transonic 
and  supersonic  efficiency  dictates  thin  airfoils  on  the 
order  of  5-8%  chord.  A  typical  military  airfoil  config¬ 
uration  along  with  nomenclature  and  the  flow  physics 
involved  is  shown  in  Fig  1. 


Fig.l  Military  high-lift  multi-element  airfoil 


Modem  airfoil  shapes  for  high  performance  aircraft 
will  need  to  be  thin  to  operate  in  the  transonic  and  su¬ 
personic  flow  regimes  efficiently.  New  integrated  de¬ 
sign  constraints  other  than  aerodynamic  performance 
criteria,  however,  will  force  design  compromises  that 
may  dictate  very  sharp  leading  edges  and  unique  plan- 


form  shapes  in  the  future.  Flap  gap  and  overhang  and 
high  lift  system  design  conducted  at  sub-scale  Rey¬ 
nolds  numbers  may  not  remain  optimal  at  flight  Rey¬ 
nolds  numbers.1  These  effects  can  become 
pronounced  in  three  dimensional  flows  due  to  the 
mechanisms  of  attachment  line  transition,  re-laminar- 
ization,  and  vortical  flow  interactions  on  military  air¬ 
craft.  Typical  aircraft  development  programs  rely  on 
sub-scale  Reynolds  number  wind  tunnel  data  to  design 
high-lift  systems,  with  empirical  corrections  for  the 
Reynolds  number  effects.  These  methods  do  not  ac¬ 
count  for  three-dimensional  effects  and  are  based  on 
conventional  designs,  and  therefore  are  not  accurate 
enough  to  guarantee  low-speed  weapon  system  perfor¬ 
mance.  There  is  motivation,  therefore,  to  develop  CFD 
methods  to  predict  maximum  lift  as  well  as  post  stall 
aerodynamic  characteristics  for  high  lift  configura¬ 
tions  accurately.  Once  validated  for  a  class  of  configu¬ 
rations  and  flow  conditions,  the  vision  is  that  CFD 
could  be  used  to  1)  complement  wind  tunnel  data,  2) 
provide  the  ability  to  optimize  component  rigging  and 
shape  at  flight  Reynolds  numbers,  3)  provide  the  abil¬ 
ity  to  correct  small-to-moderate  scale  wind  tunnel  data 
for  Reynolds  number  and  wind  tunnel  wall  and  sup¬ 
port  system  interference  effects  and  4)  provide  support 
in  designing  subcomponents. 

A  series  of  three  wind  tunnel  tests  (T378,  T385,  and 
T396)  were  conducted  in  the  NASA  langley  Research 
Center  (LaRC)  Low  Turbulence  Pressure  Tunnel 
(LTPT)  as  part  of  a  cooperative  effort  between  the  Na¬ 
vy,  McDonnel  Douglas  Aerospace  (now  Boeing)  and 
NASA  LaRC.  Surface  pressures,  forces  and  moments, 
were  measured  on  a  two  dimensional  (2-D)  airfoil 
model  of  a  Boeing  advanced  fighter  wing  section  con¬ 
figured  with  a  deflected  leading  edge  flap,  shroud  and 
a  slotted  trailing  edge  flap.  Results  from  the  previous 
tests  T378  and  T385  have  been  reported  in  the  paper 
(Ref.  2).  The  objective  of  the  present  paper  is  to  focus 
on  results  of  the  test  T396  and  Navier-Stokes  compu¬ 
tations  conducted  to  further  understand  high-lift  aero¬ 
dynamics  about  the  high-lift  fighter  airfoil  at  flight 
Reynolds  numbers  and  lift  coefficients  up  to  clmax.  Ef¬ 
fects  of  Reynolds  number,  trailing  edge  flap  gap  and 
overhang,  Gurney  flap,  vortex  generators  were  investi¬ 
gated  in  the  experiment  T396.  Some  experimental  data 


involving  transition  detection  using  hot  films  are  also 
analyzed  from  the  test  T385  and  presented  in  this  pa¬ 
per. 

Viscous  computations  were  performed  using  struc¬ 
tured/chimera  grid  for  several  of  these  configurations. 
An  incompressible  Navier-Stokes  method  was  used 
with  the  Baldwin-Barth3  (BB),  Spalart-Allmaras4 
(SA)  one  equation  turbulence  models  and  the  Shear 
stress  transport5  (SST)  two-equation  turbulence  model 
to  predict  the  Reynolds  number  effects.  Trends  in 
maximum  lift  as  a  function  of  Reynolds  number  and 
Mach  number  have  been  accurately  predicted  indicat¬ 
ing  that  the  CFD  method  could  be  used  to  extrapolate 
sub-scale  high-lift  system  aerodynamic  performance 
to  flight  scale  at  angles  of  attack  up  to  and  including 
maximum  lift6’7.  The  effects  of  flow  control  mecha¬ 
nisms  such  as  Gurney  flaps  were  simulated  using 
Navier-Stokes  analysis  and  are  compared  with  the  ex¬ 
perimental  results. 

Wind  tunnel  test  T396 


A  detailed  account  of  the  LTPT  test  facility  and  the  de¬ 
scription  of  the  high-lift  airfoil  model  is  given  in  ref. 
2.  The  model  had  a  3  ft.  span  and  22  inch  chord  and 
was  mounted  6  inches  above  the  centerline  of  the  3  ft. 
wide,  7  ft.  tall  test  section.  In  the  tests  T378  and  T385 
data  were  collected  with  the  leading  edge  flap  (LEF), 
trailing  edge  flap  (TEF)  and  the  shroud  set  at  approach 
and  clmax  conditions.  In  the  test  T396,  the  deflections 
of  the  LEF,  TEF  and  the  shroud  were  set  at  34°,  35°, 
and  22.94°  respectively  representing  a  landing  config¬ 
uration.  Pressure  and  force  data  due  to  variations  in 
gap  and  overhang  parameters  (see  table  1)  were  ob¬ 
tained.  The  baseline  configuration  was  set  with  over¬ 
hang  and  gap  of  2.66%C  and  0.512%C  respectively. 
Pressures  and  force  and  moment  data  were  obtained. 


Table  1  Parametric  overhang  and  gap  settings 


O.H.  in%C 

Gap  in%C 

2.66* 

0.512* 

1.756 

0.522 

2.756 

0.522 

Baseline 


Table  1  Parametric  overhang  and  gap  settings 


O.H.  in%C 

Gap  in%C 

3.756 

0.522 

1.756 

1.085 

2.756 

1.085 

3.756 

1.085 

1.885 

0.799 

2.885 

0.799 

3.885 

0.799 

Drag  data,  from  wake  measurements,  were  obtained 
for  selected  configurations  at  the  conclusion  of  the 
test.  In  addition  to  the  gap  and  overhang  data,  Gurney 
flap  and  vortex  generators  were  also  attached  sepa¬ 
rately  at  various  positions  on  the  high  lift  airfoil  and 
data  with  their  presence  were  obtained. 

Two  Gurney  flap  configurations  shown  in  the  figure  1. 
were  attached  on  the  lower  moldline  of  the  trailing 
edge  flap.  The  small  Gumey  flap  (GF3)  was  located 
along  the  trailing  edge  of  the  TEF.  The  large  Gumey 
flap  was  located  along  the  trailing  edge  of  the  TEF 
(GF1),  one  “pad  length”  forward  of  the  trailing  edge 
(GF4)  and  two  “pad  lengths”  forward  of  the  trailing 
edge  (GF2). 

Three  different  sized  vortex  generators  (VGs)  were 
tested  at  five  different  chord  locations.  The  0.03-in. 
high  VG1  was  located  on  the  lower  moldline  of  the 


Large  Gumey 


Small  Gumey 


_ i 

0.6% 

i 

=m]o.3% 


Fig,  1  Sketch  of  the  two  Gumey  flaps  used  in 
experiments 
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LEF  slightly  aft  (0.05-in.)  of  the  apex.  The  0.04-in. 
high  VG2  was  located  aft  of  the  shroud  knee  at  12% 
X/C.  The  0. 10-in.  high  VG3  was  located  further  aft  of 
the  shroud  “knee”  at  16.84%  of  X/C.  The  VG4 
configuration  had  two  sets  of  0.10-in.  high  VGs 
located  at  16.84%  of  X/C  and  68.98%  of  X/C.  The 
0. 10-in.  high  VG5  was  located  on  the  crown  of  the 
TEF  at  77.06%  of  X/C. 

Test  Conditions  and  Flow  Quality 

The  tests  were  conducted  at  Reynolds  numbers  from 
5, 9  and  16  million  at  the  Mach  number  of  0.2  and  at  a 
range  of  angles  of  attack.  The  large  contraction  ratio 
(17.6: 1)  of  the  LTPT  and  nine  anti-turbulence  screens 
provide  excellent  free  stream  flow  quality.  During  the 
NASA  LTPT  entries,  Test  378,  Test  385  and  Test  396, 
it  was  important  to  develop  2D  flow  inside  the  tunnel. 
Side  wall  boundary  layer  control  system  helped  to 
ensure  2-D  flow.  The  real  time  observation  of 
sufficiently  “flat”  traces  of  the  pressures  measured  by 
the  spanwise  taps  assured  2-D  flow.  When  the 
spanwise  pressures  deviate  beyond  the  acceptable 
tolerance  level,  the  side  wall  suction  was  adjusted  to 
ensure  2-D  flow,  side  wall  venting  equal  to 
approximately  0.2%  of  the  freestream  mass  flow  was 
used  to  remove  side  wall  boundary  layers.  Since  the 
side  wall  vents  are  metric,  there  was  a  lack  of 
confidence  in  the  accuracy  of  the  measured  balance 
data.  Therefore,  to  obtain  lift,  the  pressures  measured 
by  surface  pressure  taps  was  integrated.  Lift  was  not 
corrected  for  wind  tunnel  effects  since  a  previous 
study  investigating  wall  effects  on  a  commercial 
configuration  concluded  that  these  effects  were 
minimal  at  lift  coefficients  below  2.5  . 

Results  of  the  Wind  Tunnel  Experiments 

Effects  of  Gap  and  overhang.  Gumev  flaps  and  Vortex 
generators 

The  basic  results  of  the  experiments  are  presented  in 
this  abstract  and  the  complete  results  will  be  dis¬ 
cussed  in  detail  in  the  actual  paper.  The  effects  of  the 
variation  in  the  gap  and  overhang,  vortex  generators 
and  Gumey  flap  at  the  Reynolds  number  of  5  and  16 


million  were  documented  in  detail.  Out  of  several 
overhang  and  gaps,  overhang  and  gap  of  1.085%  and 
1.756%  respectively  provided  the  maximum  increase 
of  3%  in  clmax  over  the  baseline  case  as  seen  in  the 
figure  2.  Attaching  vortex  generators  on  the  flap  pro¬ 
duced  slightly  lower  benefits  than  the  Gurney  flap  or 
the  optimum  overhang  and  gap  arrangements.  The 
maximum  increase  in  clraax  due  to  the  large  Gumey 
flap  on  the  lower  surface  of  the  TEF  is  shown  also 
shown  in  the  same  figure.  The  Gumey  flaps  were 
attached  to  the  Baseline  configuration.  The  Gumey 
flaps  produced  increase  in  clmax  of  about  4%  which 
amounts  to  roughly  about  2%  reduction  in  approach 
speed. 


Fig.  2  The  effects  of  gap  and  overhang,  Gumey 
flap  and  vortex  generators  on  q 

Hot  Film  Measurements  • 

In  order  to  understand  the  flow  physics,  hot  film  gaug¬ 
es  were  attached  near  the  leading  edge  regions  of  the 
wing  and  the  TEF,  and  upper  surface  of  the  wing.  Four 
different  hot  film  patches  were  used  with  slightly  dif¬ 
ferent  reference  conditions.  RMS  signals  were  mea¬ 
sured  and  correlation  coefficients  were  determined 
using  hot  film  measurements  for  a  range  of  angle  of  at¬ 
tack  and  for  several  Reynolds  numbers  and  flap  defec¬ 
tions. 


Navier-Stokes  Solutions 


The  RMS  values  and  the  correlation  coefficients  are 
presented  in  Fig.  13  for  a  Reynolds  number  of  16  mil¬ 
lion  and  at  AO  A  of  1°.  The  LEF,  shroud  and  TEF  are 
deflected  at  34°,  35°  and  22.94°  respectively.  Typical¬ 
ly,  RMS  values  will  be  low  in  the  laminar  boundary 
layers,  reach  a  peak  at  transition,  and  then  will  de¬ 
crease  to  a  value  lower  than  peak  transition  but  sub¬ 
stantially  higher  than  that  of  a  laminar  boundary  layer. 
The  leading  edge  region  seems  to  be  characterized  by 
a  very  short  bubble  at  the  leading  edge  followed  by 
transition  due  to  laminar  separation.  However,  it  is  to 
be  noted  that  in  the  computational  analysis  the  entire 
flow  was  assumed  to  be  turbulent.  A  detailed  analysis 
of  the  hot  film  measurements  and  their  correlation 
with  CFD  data  and  flow  physics  will  be  given  in  the 
actual  paper. 


Fig.  3  RMS  Signals  and  Correlation  coefficients  from  hot 
film  measurements  for  AOA=l°,  Reynolds  number  of  16 
million. 


Viscous  computations  were  performed  using 
structured  /  chimera  grid  for  several  of  these 
configurations  with  the  Baldwin-Barth,  Spalart- 
Allmaras  one  equation  turbulence  models  and  the 
Shear  stress  transport  two-equation  turbulence  model. 


Grid  generation  and  Solution  Method 

The  Chimera  or  overset  grid  used  for  the  computations 
is  shown  in  Fig  4.The  field  grid  was  generated  by  us- 


Fig.  4  Chimera  multi-element  airfoil  initial  grid 


ing  a  hyperbolic  grid  generator  for  the  main  element 
field  grid  and  algebraic  and  elliptic  grid  generators  for 
all  other  components.  The  first  layer  off  the  surface 
was  specified  at  a  normal  distance  of  1 .0  x  10‘6  chords 
for  each  element.  This  spacing  resulted  in  y+  values 
equal  to  or  less  than  1  for  a  freestream  Reynolds  num¬ 
ber  of  16  million.  A  total  of  31  points  were  clustered 
within  the  boundary  layer  over  most  of  the  configura- 
tion.The  total  number  of  points  in  the  initial  grid  was 
71,630.  The  main  element  grid  extended  to  the  far- 
field  which  was  placed  at  a  minimum  distance  of  25 
chords  away  from  the  configuration  in  all  directions. 
The  wakes  of  the  main  element  and  flap  were  expand¬ 
ed  downstream  to  reduce  aspect  ratio  of  cells  at  wake 
computational  planes.  Holes  were  cut  and  connectivi¬ 
ty  using  single  fringes  between  grids  was  specified  by 
using  the  GMAN  module  of  MDA’s  grid  system  soft¬ 
ware 

The  ENS2D9  incompressible  Navier-Stokes  flow  solv¬ 
er  was  used  in  the  analysis.The  codes  have  been  used 
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extensively  by  the  commercial  industry  to  predict 
high-lift  airfoil  flows.  The  code  INS2D  uses  the  artifi¬ 
cial  compressibility  approach  to  couple  the  mass  and 
momentum  equations  which  generates  a  hyperbolic 
set  of  equations.  It  employs  a  third-order  accurate  up¬ 
wind  formulation  for  the  resulting  convective  terms.  A 
recent  enhancement  to  the  solver  is  the  incorporation 
of  a  GMRES  iterative  algorithm  which  reduces  ran 
time  significantly.10  Implemented  turbulence  models 
include  the  BB  and  SA  one-equation,  models  and  the 
SST  two-equation  model.  The  empirical  parameters 
used  in  the  SA  model  are  adjusted  slightly  as  de¬ 
scribed  in  Rogers1 1  The  code  is  capable  of  reading  in  a 
2D  version  of  PEGASUS  style  grid  inter-connectivity 
information. 

Grid  Refinements 

Since  there  were  differences  in  boundary  conditions 
used  and  substantial  gradients  in  areas  of  flow  over 
coarse  grid  points  such  as  at  the  leading  edge,  the  grid 
was  refined.  To  reduce  the  effects  of  far-field  boundary 
conditions,  a  Cartesian  background  grid  was  con¬ 
structed  extending  50  chord  lengths  away  from  the  air¬ 
foil  in  all  directions.  Grid  points  were  clustered  in  the 
streamwise  directions  at  the  leading  edge,  the  leading 
edge  knee,  the  shroud  hinge  line,  the  trailing  edges  of 
the  main  element  and  trailing  edge  flap  element,  and  in 
the  ledge,  cove,  and  slot  flow  areas.  Grid  clustering 
normal  to  the  walls  remained  unchanged.  An  example 
of  the  dense  grid  clustering  at  the  leading  edge  and 
shroud  trailing  edge  is  shown  in  Fig.  5. 


clustering  of  refined  grid. 


The  location  of  the  computational  wake  boundary  aft 


of  the  main  element  was  moved  to  align  with  the  cen¬ 
ter  of  the  wake  at  a  =  0°.  The  total  number  of  grid 
points  in  the  refined  grid  was  85,957.  Although  there 
were  substantial  improvements  in  grid  clustering,  the 
refined  grid  had  only  a  modest  increase  in  total  grid 
points. 

Results  of  the  Computational  Analysis 

The  flow  over  both  initial  and  refined  grids  was  solved 
using  the  INS2D  code  at  a  Reynolds  number  of  16  mil¬ 
lion,  at  angles-of-attack  between  a  =  -10°  and  a  =  1 1°, 
using  the  BB,  SA,  and  SST  turbulence  models.This 
Reynolds  number  corresponds  to  actual  flight  condi¬ 
tions  for  a  section  of  a  modem  Navy  strike/fighter  air¬ 
craft  wing  during  take-off  and  landing.  All 
computations  imposed  fully  turbulent  flow  every¬ 
where,  and  convergence  was  considered  to  be  five  sig¬ 
nificant  digits  of  the  lift  coefficient. 

Assesment  of  turbulence  models 

In  figure  6,  lift  coefficients  are  presented  using  the  BB, 
SA,  and  SST  models  on  the  refined  grid.  The  resulting 
lift  coefficient  for  all  the  three  turbulence  models  were 
nearly  identical  at  flow  conditions  below  clmax.  For 
post  stall  lift,  the  BB  and  SST  models  predicted  nearly 
the  same  and  had  the  same  trends  as  the  experiment. 
The  SA  model  did  not  predict  the  drops  in  experimen¬ 
tal  lift  at  a  =  6°,  and  a  =  1 1°.  The  difference  in  comput¬ 
ed  lift  at  an  AOA  of  a  =  -6°  is  caused  by  the  differences 
in  separation  and  shear  layer  trajectory  predicted  by 
the  BB  and  SA  models,  respectively,  as  the  shear  layer 
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convects  below  the  airfoil. 


Fig.  6  Lift  coefficcient  for  refoned  grid  using  BB. 
SA  and  SST  turbulence  models 


Effect  of  Gumev  Flaps 

Computational  analysis  were  done  with  a  2%  Gumey 
flap  on  the  lower  surface  at  the  trailing  edge  of  the  TEF 
with  Baldwin-Barth  turbulence  model  and  are  shown 
in  Fig.  7.  Higher  suction  was  obtained  on  the  upper 
surface  due  to  the  Gumey  flap  in  the  computations.  As 
seen  in  the  figure  right  trend  in  Cj  is  predicted  with  the 
computations. 


Fig.  7  Lift  coefficient  with  2%  Gumey  flap  with  BB 
turbulence  model  at  16  million  Reynolds  number 


Acknowledgments 

The  work  was  sponsored  by  the  ONR  Air  Vehicle 
Technology  Program.  Some  of  the  computations  were 
conducted  at  the  NAS  facility  at  NASA  Ames  and  at 
the  CEWES  HPC  MSRC  in  Mississippi. 

References 

1.  Meredith,  P.  T.,  “Viscous  Phenomena  Affecting 
High-Lift  Systems  and  Suggestions  for  Future  CFD 
Development,”  AGARD-CP-515,  September  1993, 
pages  19-1,  19-8. 

2.  Hobbs,  C.R,  Spaid,  F.,  Ely,  W.,  and  Goodman,  W., 
“High  Lift  Research  Program  for  a  Fighter-type, 
Multi-element  Airfoil  at  High  Reynolds  Numbers,” 
AIAA  96-0410,  January  1996. 

3.  Baldwin,  B.  S.,  and  Barth  T.  J.,  “A  One-Equation 
Turbulence  Transport  Model  for  High  Reynolds 
Number  Wall  Bounded  Hows”  NASA  TM  102847, 
August  1990. 

4.  Spalart,  PR.  and  Allmaras,  S.R.,  “A  One-Equation 
Turbulence  model  for  Aerodynamic  Hows”  AIAA 
92-0439,  January  1992. 

5.  Menter,  F.R.,  “Zonal  Two  Equation  k-w  Turbulence 
Models  for  Aerodynamic  Hows,”  AIAA  93-2906, 
July  1993. 

6.  Kem,  S.  B.,  “Evaluation  of  Turbulence  Models  for 
High-Lift  Military  Airfoil  Flowfields”,  AIAA-96- 
0057,  January  1996. 

7.  Hall,  D.,  and  Kem,  S.,  “  Prediction  of  REynolds 
number  Effects  of  Military  High-Lift  Airfoil”,  AIAA 
96-2400,  June  1996. 

8.  Cao,  H.V.,  Kusunose,  K.,  Spalart,  P.R.,  Ishimitsu, 
K.K.,  Rogers,  S.E.,  and  McGhee,  R.J.,  “Study  of 
Wind  Tunnel  Wall  Interference  for  Multi-Element 
Airfoils  Using  a  Navier-Stokes  Code,”  AIAA  94- 
1933,  June  1994. 

9.  Rogers,  S.E.,  and  Kwak,  D.,  “An  Upwind  Differ¬ 
encing  Scheme  for  the  Incompressible  Navier-Stokes 
Equations,”  NASA  TM  101051,  November  1988,  also 
appeared  in  Applied  Numerical  Mathematics,  Vol.  8, 
1991,  pp.  43-65. 

10.  Rogers,  S.E.,  “A  Comparison  of  Implicit  Schemes 
for  the  Incompressible  Navier-Stokes  Equations  with 
Artificial  Compressibility,”  AIAA  95-0567,  January 


1995. 

11.  Rogers,  S.E.,  Menter,  F.R.,  Durbin,  P.A.  and  Man- 
sour,  N.N.,  “  A  comparison  of  Turbulence  Models  in 
Computing  Multi-Element  Airfoil  Flows, “  AIAA  94- 
0291,  January  1994. 


